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*vj Abstract 

t-H It is well known that a free three-dimensional rigid body admits three stationary rota- 

tions. These are the rotations around three principal axes of inertia. The rotations around 
the long and the short axes are stable, while the rotation around the intermediate axis is 
unstable. We generalize this result to the case of a multidimensional body. The stability 
problem is being solved for a dense subset of stationary rotations in arbitrary dimension. 
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1 Introduction 



The above treatment of the problem of 
rotation may, in contradistinction to the 
usual method, be transposed, word for 
word, from three-dimensional space to 
multidimensional spaces. This is, indeed, 
irrelevant in practice. On the other 
hand, the fact that we have freed 
ourselves from the limitation to a 
definite dimensional number and that we 
have formulated physical laws in such a 
way that the dimensional number 
appears accidental in them, gives us an 
assurance that we have succeeded fully in 
grasping them mathematically. 

H.Wcyl. Space, Time, Matter. 



1.1 Statement of the problem 

Speaking informally, free multidimensional rigid body is simply a rigid body rotating in 
multidimensional space without action of any external forces (i.e. by inertia). 

Let us first discuss a three-dimensional free rigid body (the so-called Euler case in the 
rigid body dynamics). A good model for such a body is a book or a parallelepiped shaped 
box. 

Throw a book in the air spinning it in arbitrary direction. If we neglect the gravity force, 
then what we get is exactly the Euler case. Note that a general trajectory of the book is not 
a rotation in the usual sense. At each moment of time the book is indeed rotating around 
some axis, but this axis is changing as time goes. What we are interested in, are the relative 
equilibria of the system, i.e. such trajectories for which the axis of rotation remains fixed. 
Such rotations are also called stationary. 

It is well known that a rotation of a generic three-dimensional rigid body (i.e. a body 
with pairwise distinct principal moments of inertia) is stationary if and only if the axis of 
the rotation is one of the three principal axes of inertia. If we deal with a homogeneous 
parallelepiped shaped body, then the principal axes of inertia coincide with the axes of 
symmetry. 

Let S > I > L be the principal moments of inertia. The corresponding principal axes of 
inertia are called, respectively, the short, the intermediate, and the long axis of inertia. In 
the case of a symmetric body these names correspond to the actual lengths of the axes of 
symmetry (see Figure [I]) . 




Figure 1: Three-dimensional rigid body. 



Stationary rotations around different principal axes of inertia have different dynamical 
features (see [IH2])- To demonstrate it, spin the body around one of these axes. Of course, 
since our hands and eyes are not too precise, it is not possible to spin exactly around the 
chosen axis. But this is not going to be fatal if this axis is the short or the long one. We will 
not have a stationary rotation, but something very close to it. This is due to the fact that the 
rotations of a free three-dimensional rigid body around the short and the long principal axes 
of inertia are (Lyapunov) stable. But if we spin around the intermediate axis, we will see 
something essentially different: the axis of rotation will start changing rapidly and the body 
will start rotating in other direction. This is because the rotation of a free three-dimensional 
rigid body around the intermediate principal axis of inertia is (Lyapunov) unstable. 

Basically, what we are interested in, is to generalize this result to the case of a multidi- 
mensional body. 

Remark 1.1. The system describing the motion of a free multidimensional rigid body (see 



Section 1.3) is also known as the Manakov top. Manakov's achievement was to prove that 
this system is integrable (see [4])- However, the equations themselves go back to Frahm [5]. 
fn the famous paper [B] Arnold wrote these equations in the form of hamiltonian equations 
on so(n)* and generalized them to the case of an arbitrary Lie algebra (see also [2]). A 
possibility to generalize the equations of a free rigid body to the multidimensional case was 
also mentioned by Weyl [7J. 

First we shall discuss how an n-dimensional body may rotate. At each moment of time 
K n is decomposed into a sum of m pairwise orthogonal two-dimensional planes fli, . . . , H m 
and an n — 2m-dimensional space flo orthogonal to all these planes: 



= (© n <) ® n °- 



There is an independent rotation in each of the planes ITi, . . . ,n m , while Ilo is fixed. This 
is just a reformulation of the theorem about canonical form of a skew-symmetric operator. 
Note that flo may be zero in the even-dimensional case, which means that there are no fixed 
axes. 

A rotation is stationary if all the planes Flo, ... , IT m don't change with time (this condi- 
tion automatically implies that the velocities of the rotations are also constant). 

Before studying stationary rotations for stability, we need to find these rotations. Recall 
that a rotation of a generic three-dimensional rigid body is stationary if and only if it is 
a rotation around one of the principal axes of inertia. In the multidimensional case the 
situation is slightly more complicated. If the planes Ilo, . . . ,n m are spanned by principal 
axes of inertia (such rotations are called in [8] regular), then the rotation is stationary. But 



the converse is not necessarily true (see [5] and Section 1.5 of the present paper). However, 
it is true provided that the angular velocities of the rotations in the planes Hi , . . . Xl m are 
pairwise distinct (see [9] HQ]). Moreover, we will show that the rotations for which the 
planes Ho, . . . ,n m are not spanned by principal axes of inertia (such rotations are called 
in [8] exotic) are always unstable. Therefore, we will mainly deal with regular stationary 
rotations. 

In four dimensions the problem of stability was studied from different points of view by 
Oshemkov 11 , Feher and Marshall 9 , Birtea, Ca§u, Ratiu, and Turhan [10] . Birtea and 
Casu [12] . The answer is known for a dense subset of relative equilibria. 

Remark 1.2. Oshemkov didn't actually pose the stability problem. However, the answer 
can be easily obtained from his bifurcation diagrams. 

The five-dimensional case was studied by Ca§u in [13] • The set of equilibria which were 
studied in this case is not dense. 

General even-dimensional rigid body was discussed in Spiegler's PhD thesis [14]. A 
sufficient condition for an equilibrium to be stable is found. 



In this paper we solve the stability problem for a dense subset of relative equilibria in 
arbitrary dimension. This result is presented in Section [2] 

Recall that in three dimensions stability depends only on the rotation axis. In many 
dimensions the situation is slightly more complicated: angular velocities of rotations may 
also affect stability. However, the answer is still rather elegant, and the classical three- 
dimensional result may be viewed as a particular case of this general answer. 

1.2 Methods of stability investigation 

In this section we will discuss different approaches to the problem of stability of relative 
equilibria of a multidimensional rigid body. 

Spiegler [14] approached the problem using the method known as Arnold or energy- 
Casimir method (the method was introduced in [6|; see also [2]). This method can be 
formulated as follows: Let x be an equilibrium point of a hamiltonian system on a Poisson 
manifold. Then a; is a critical point for the restriction of the Hamiltonian to the symplectic 
leaf O passing through x. If this critical point is a non-degenerate minimum or maximum, 
and O is regular, then x is stable. 

The energy-Casimir method is a very powerful tool for studying equilibria of general 
Hamiltonian systems (especially inifinite-dimensional, see [3j I15H17] for an overview) . But 
if a system possesses some additional symmetries (for example, it is completely integrable), 
there are more efficient ways to study stability. Let us consider the following simple example: 
take the Hamiltonian 

n 

If there are cj^'s of different sign (which is very likely if n is large), then the origin is not 
a minimum or a maximum point for H. Consequently, we can't apply the energy-Casimir 
method. But if we observe that the system possesses an integral 

n 

/ = ^pl + qh 

then stability is immediately established. 

This example shows what normally happens in integrable systems. Since the Birkhoff 
normal form for an integrable system always converges (see |18l I19| ). the Hamiltonian can 
be brought to the form M in the neighborhood of a generic linearly stable equilibrium. 
Since cj^'s are, in general, of different signs, the energy-Casimir method in the integrable 
case does not work for "most" stable equilibria. 

As it was shown by Manakov in [4], the system describing the motion of a multidimen- 
sional rigid body is integrable. Therefore we may apply methods developed in the theory of 
integrable hamiltonian systems. 

Remark 1.3. More precisely, Manakov showed that the system admits an L — A pair with 
a spectral parameter, which allowed him to find integrals. However, he did not prove these 
integrals are enough for Liouville integrability. Complete integrability in the Liouville sense 
was proved by Fomenko and Mischenko in 20 and by Ratiu in |21j . See also earlier paper 
[22] . where the quadratic integrals of the problem are found, and [23], where these integrals 
are proved to be in involution. 

Manakov also showed that the system is integrable in ^-functions of Riemann surfaces 
(see also [53], where the four-dimensional case is integrated by elementary means). 

Let us mention two possible ways to study stability of equilibria of an integrable hamil- 
tonian system: 

1) Extended energy-Casimir method. The idea is to replace the Hamiltonian with any 
other integral in the formulation of the energy-Casimir method. The extended method can be 



formulated as follows: Let x be an equilibrium point of a hamiltonian system on a Poisson 
manifold, and let O be the synrplectic leaf passing through x. Suppose that the system 
possesses integrals fi,...f n . Assume that we can find a linear combination / = X} <-*/*> 
such that the restriction of / to O has a non-degenerate minimum at the point x. Then, if 
O is regular, x is stable. 

Remark 1.4. Note that if x is not critical for the restriction of / to O, then the whole 
trajectory of the hamiltonian flow generated by / passing through x consists of equilibrium 
points. In this situation it is possible to show that x is unstable, provided the system is 
non-resonant. 



This method was applied in four dimensions in 9, 12 . In five dimensions it was applied 
in [13], however only quadratic integrals were considered (while this system possesses also a 
cubic integral). 

2) Eliasson theorem (this theorem is actually due to Russmann [25], Vey [26], Eliasson 
27, 28 and Miranda 29 ). The theorem provides a normal form for the momentum mapping 
in the neighborhood of any equilibrium satisfying the so-called non-degeneracy condition. 
There is only a finite number of normal forms. They are distinguished by the so-called 
Williamson type of a point. Only one type corresponds to stable behavior. Therefore, the 
stability problem is easily solved provided we can verify non-degeneracy and determine the 
type of a point. See [30] for precise defintions. 

The disadvantage of this method is that it works only for non-degenerate points. Pro- 
vided the non-degeneracy condition is satisfied, the method is equivalent to the extended 
energy-Casimir method. However, if it is not satisfied, it is still possible that the energy- 
Casimir method works, while the Eliasson theorem does not. 

The advantage of the method is that we get much more information, than simply stabil- 
ity/instability. 

This method was applied in four dimensions in [10] . However, all the necessary informa- 
tion can be extracted from the bifurcation diagrams constructed in [11| . 

Both methods are very efficient, but their explicit application involves heavy calculations 
if the dimension is larger then four, because the number and the degree of integrals grow 
with the dimension. This brings implicit methods to the stage. One of such methods is to 
make use of the presence of a bihamiltonian structure. 

Remark 1.5. The notion of a bihamiltonian system was introduced by F. Magri in |31j . 
Bihamiltonian structure is known to be closely related with integrability. On the one hand, 
if a system is bihamiltonian, then it admits "many" integrals in involution (see [31-34 ). On 
the other hand, many integrable equations coming from mechanics, mathematical physics 
and geometry possess a bihanriltonian structure (see, for example, [311 1354443] ). 

The bihamiltonian structure for the multidimensional rigid body equations was discov- 
ered by A.Bolsinov in [44] (see also [45H47| ) . This structure is defined on the dual space of 
the Lie algebra of skew-symmetric matrices. 

Remark 1.6. However, it is possible to give another bihamiltonian formulation: the bi- 
hamiltonian structure is defined on the dual of sl(n)* and then the system is obtained by the 
restriction from sl(n)* to so(n)* . This structure is implicitly present in [20] and explicitly 
written in gBj. See also [21] . 

The bihamiltonian approach for studying relative equilibria of a multidimensional rigid 
body was introduced in [35]. In this paper (among other results) Bolsinov and Oshemkov 
obtain a sufficient condition for non-degeneracy. Developing their ideas and applying the 
technique of [50], we give a necessary and sufficient condition for an equilibrium to be non- 
degenerate and determine type of all non-degenerate points. This allows us to solve the 
stability problem for a dense subset of relative equilibria. 



1.3 The equations 

The motion of a free multidimensional rigid body is described by the Euler-Arnold (or 
Euler-Frahm) equations on so(n)* (identified with so(n)). These equations have the form 

M =[M,Q] 

where 

• M € so(n)* is a skew-symmetric matrix, called the angular momentum matrix. 



• J is symmetric. We call it the covariance matrix of the body (see Section 1.4 1. 

• H is a skew-symmetric matrix, called the angular velocity matrix. It is uniquely defined 
by the relation 

m = nj + Jfi. 

For explicit derivation of these equations see |21j . 

Remark 1.7. Since the map J: so(n) — > so(n) given by the formula 

j(n) = nj + jn (3) 

is invertible, our equations can be rewritten in the ^-coordinates: 

ti = j- 1 ([j(n),a\). 

However, the explicit formula for J~ x is complicated, therefore it is convenient to introduce 
the variable M and write down the equations in the form H. 

Remark 1.8. Note that the equations (pi describe only the dynamics of the angular velocity 
matrix. If we want to recover the dynamics in the whole phase space T*SO(n), we should 
add Poisson equations 

x = xn. 

However, we will only be interested in the reduced dynamics, given by the system B. Note 
that relative equilibria of a rigid body are nothing else but the equilibrium points of the 
system Q2p. 

1.4 The inertia tensor versus the covariance matrix 

In the multidimensional case the matrix J entering the equations H is sometimes referred 
as the "inertia tensor". This is not very precise, because in the three-dimensional case 
the inertia tensor is not J, but the map J: so(3) — > so(3) given by the formula pi. The 
coincidence of dimensions of R 3 and so (3) makes it possible to consider both J and J as 
operators (or bilinear forms) on R 3 . In the multidimensional case these operators act on 
different spaces. We suggest the following definitions: 

• The inertia tensor is J . J is the covariance matrix (tensor) of the body. Recall (see 
|21| ) that the entries of J are the values of the covariance of the coordinates: 

Jij = (xi- Xi)(xj - Xj)dn, 

where Xi are the center of mass coordinates. 

• The principal moments of inertia are the eigenvalues of the inertia tensor J . If we 
denote the eigenvalues of the covariance matrix J by A;, then the principal moments 

Xi + Xj for i -£ j. 



• The principal axes of inertia are the eigenvectors of the covariance matrix. If we denote 
them by e,, the eigenvectors of the inertia tensor J will be d A e? for i 7^ j. Note that 
in the three-dimensional case the eigenvectors of the inertia tensor and the covariance 
matrix coincide. 

• A body is generic if the eigenvalues of the covariance matrix are pairwise distinct. If 
all the principal moments of inertia are pairwise distinct, then the body is generic. 
The converse is not necessarily true if the dimension is larger than three. In three 
dimensions these two statements are equivalent. 

1.5 Description of relative equilibria 

Here we briefly state the results of [8J. 

Theorem 1. Consider the system of Euler- Arnold equations pi). Suppose that J has pair- 
wise distinct eigenvalues. Then M is an equilibrium point of the system if and only if there 
exists an orthonormal basis such that J is diagonal, and Q. is block-diagonal of the following 
form 

UkAh 



Q = 



(4) 



V ■• / 



where Ai £ so(2mi) n SO(2m;) for some mt > 0, and uji '$ are distinct positive real numbers. 

Corollary 1.1. Suppose that M is a relative equilibrium, all eigenvalues of J are pairwise 
distinct. Moreover, let all non-zero eigenvalues of Q (eigenfrequencies of rotation) be pair- 
wise distinct. Then there exists an orthonormal basis such that J is diagonal, while Q and 
M are block-diagonal with two-by-two blocks on the diagonal. 

In other words, a stationary rotation with pairwise distinct eigenfrequencies is a rotation 
"in principal axes of inertia". 

Definition 1. We will say that an equilibrium M is regular if there exists an orthonormal 
basis such that J is diagonal and Q is block-diagonal with two-by-two blocks on the diagonal 
(i.e. M is a rotation "in principal axes of inertia"). Otherwise, we will say that M is exotic. 

Corollary |1.1| says that all stationary rotations with pairwise distinct eigenfrequencies 
are regular. 

We will see later that regular equilibria are exactly rank zero singular points of the 
system in the sense of integrable systems theory, which means that these points are critical 
points for all integrals. Exotic equilibria are, on the contrary, non-zero rank singular points, 
which means that there exists an integral the differential of which doesn't vanish in the 
equilibrium point. This implies that exotic equilibria are not isolated on a coadjoint orbit, 
but form whole smooth submanifolds of equilibrium points, while regular equilibria are 
always isolated (on a coadjoint orbit). 

Remark 1.9. However, it is possible for a regular equilibrium to belong to the closure of 
the set of exotic equilibria. 

We will also show that exotic equilibria are always unstable. 



2 Stability 

2.1 Parabolic diagram of a regular relative equilibrium 

Suppose we have a regular relative equilibrium. Then there exists an orthonormal basis such 
that J is diagonal and fl is block-diagonal with two-by-two blocks on the diagonal. In other 
words, there exists a decomposition 

/ m \ 

011; ©n , (5) 



where 

• Each II;, i > is spanned by two principal axes of inertia. There is a rotation in each 
of these planes with angular velocity w; . 

• ITo is spanned by n — 2m principal axes of inertia and is fixed. 

Therefore, to define a regular relative equilibrium, we need to choose an integer m such that 
< m < [n/2] and pick m pairs out of the set of the principal axes of inertia. For each 
pair we need to define an angular velocity. Choice of pairs and values of angular velocities 
uniquely define a regular relative equilibrium. Knowing this data, we want to understand 
whether an equilibrium is stable or not. 

We are going to define an object called the parabolic diagram of a regular relative equi- 
librium. This object will allow us to express stability conditions in geometric terms. 

For each plane II;, i > let us denote by Ai (II;), A2(fl;) the eigenvalues of the covariance 
matrix J corresponding to the principal axes of inertia which span IT;. By w(IIi) denote the 
angular velocity of rotation in the plane II;. 

Definition 2. The parabolic diagram of a regular relative equilibrium is the following set 
of parabolas and vertical lines drawn on the same coordinate plane: 

• For each EI; draw a parabola given by 

v M = (x - X 1 (U i ) 2 )(x - \2JUj) 2 ) 

XA ' a;(no 2 (Ai(no + A 2 (nO) 2 • { °> 

• For all fixed principal axes draw vertical lines through the squares of corresponding 
eigenvalues of J. 

Remark 2.1. Each Xi is a quadratic function the roots of which are the squares of the 
eigenvalues of J. Therefore what we do is we simply draw parabolas through the squares 
of eigenvalues corresponding to moving axes and vertical lines through the squares of eigen- 
values corresponding to fixed axes. 

Remark 2.2. The leading coefficient of Xi is inverse proportional to the square of the 
angular momentum m(IIj) = w(II;)(Ai(II;) + As (II;)). 
Let us accept the following formal agreement: 

• Two parabolas intersect at infinity, if their leading coefficients are equal, i.e. if they 
have only one point of intersection (of multiplicity one) or no points of intersection 
(neither real, nor complex). 

• Two parabolas are tangent at infinity if they have no points of intersection (real or 
complex), i.e. if they can be obtained from each other by a vertical shift. Note that 
this is only possible if the covariance matrix eigenvalues satisfy A; + A 2 = A 2 . + A 2 for 
some distinct i,j,k,l. 

After accepting this agreement, the following becomes true: 

Proposition 2.1. Any two parabolas on a parabolic diagram intersect exactly at two points 

with multiplicities. 

Definition 3. We will say that a parabolic diagram is generic if all intersections on it are 

simple, i.e. it contains no multiple intersections and no points of tangcncy. 



2.2 Stability theorems 

The following three theorems are the main results of the paper. 

Theorem 2. Let a generic multidimensional rigid body be in a state of regular relative 



equilibrium. Assume that 

• All intersections on the parabolic diagram of the equilibrium are either real and belong 
to the upper half-plane or infinite. 

• The parabolic diagram is generic. 

• The equilibrium has no more than two fixed axes. 

Then M is stable. 

Theorem 3. Let a generic multidimensional rigid body be in a state of regular relative 
equilibrium. Assume that 

• There is at least one intersection on the parabolic diagram of the equilibrium which is 
either complex or belong to the lower half-plane. 

• The equilibrium has no more than two fixed axes. 

Then M is unstable. 

Theorem 4. All exotic relative equilibria of a generic multidimensional rigid body are un- 
stable. 

Conjecture (Stability criteria). A relative equilibrium is stable if and only if it is regular 
and all intersections on its parabolic diagram are either real and belong to the upper half- 
plane or infinite. 

Whether or not this conjecture is true, Theorems [2] [3] HI are already enough to solve 
the stability problem for an open dense subset of relative equilibria. Proof of these three 
theorems can be found in Section [5] 

2.3 Examples 

Now we shall discuss some examples. First let us recover the classical three-dimensional 
result. 

Example 2.1 (Three-dimensional body). Parabolic diagrams for a 3-dimensional body are 
illustrated on Figure [2] We see from the parabolic diagrams that the rotation around the 
intermediate axis (diagram 2) is unstable, while the rotations around the short (diagram 1) 
and the long (diagram 3) axes are stable. 






1. Short axis: stable 2. Intermediate axis: unstable 3. Long axis: stable 

Figure 2: Parabolic diagrams for a 3d-body 



Example 2.2 (Four-dimensional body). Parabolic diagrams for a 4-dimensional body are 
illustrated on Figured] We assume that A? + A4 ^ Al + A|, where Ai < A2 < A3 < A4 are 
the eigenvalues of the covariance matrix. 



A four-dimensional body has three different modes of rotation: 

1. The first plane of rotation is spanned by two short axes of inertia, the second plane is 
spanned by two long axes of inertia. We look at the parabolic diagrams (see diagrams 
1,2) and see that such rotation is stable. 

2. The first plane of rotation is spanned by the shortest and the second longest axes, the 
second plane is spanned by the longest and the third shortest axes. We look at the 
parabolic diagrams (see diagrams 3,4) and see that such rotation is unstable. 

3. The first plane of rotation is spanned by the longest and the shortest axes, and the 
second plane is spanned by two intermediate axes. See diagrams 5-11. We see that 
stability depends in this case on the ratio of angular velocities. If the rotation in the 
"inner" (i.e. spanned by the intermediate axes) plane is fast enough, we have stability. 
If it is slow, we have instability. The situation is similar to the rotation of a 3d body 
with a gyroscope inside. If the gyroscope is rotating fast enough, it stabilizes the 
rotation of the body around the intermediate axis of inertia. The "inner" plane in a 
four-dimensional body plays the role of a gyroscope. 

In the non-generic case illustrated on diagram 7 our theorems do not give an answer. 
However, it is possible to show that the corresponding rotation is stable. This implies 
that stability loss in a four-dimensional body is always "soft" , i.e. the set of stable 
equilibria is closed. 




1. Stable 



2. Stable 



3. Unstable 





4. Unstable 



5. Stable 



6. Stable 




7. Unknown 



Unstable 



9. Unstable 




10. Unstable 



Figure 3: Parabolic diagrams for a 4d-body 
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3 Bihamiltonian structure and complete integrabil- 
ity 

3.1 The bihamiltonian structure 

We will denote the standard Lie bracket on so(n) by [ , ]<» an d the corresponding Lie-Poisson 
bracket on so(n)* by { , }ao. It is given by 

{f,g}^(M) = (M, [d/,d5]oo>, where M € So(n)*. 

The following is standard 

Proposition 3.1. The equations pi) are hamiltonian with respect to the Lie-Poisson bracket 
{ , }oo- The Hamiltonian is given by the kinetic energy 

H=\{£l,M) 

Now we introduce a second operation on so(n) defined by 

[X,Y] = XJ 2 Y -YJ 2 X. 
Proposition 3.2. 

1. [,}a is a Lie bracket compatible with the standard Lie bracket. In other words, any 
linear combinations of these brackets define a Lie algebra structure onso(n). 

2. The corresponding Lie-Poisson bracket { , }o on so(n)* given by 

{f,g}o = {M,{df,dg] ). 
is compatible with the Lie-Poisson bracket { , }oo- 

Consequently, we have a Lie pencil on so(n) and a Poisson pencil on so(ra)*. We will 
write down these pencils in the form 

[X, Y] x = [X, Y] - \[X, Y] x = X(J 2 - \E)Y - Y{J 2 - XE)X, (7) 

{/, 9}x = {/, 9}o -Hf,g}oo = (M, df(J 2 - \E)dg - dg(J 2 - XE)df) . (8) 

Theorem 5 (A.Bolsinov, |44H46| ). The system M) is hamiltonian with respect to any bracket 
{ ,}\, i.e. it is bihamiltonian. The hamiltonian is given by 

H A = -i((J + v / A£)" 1 0(J + \ / AS)- 1 ,M). (9) 

Remark 3.1. Since J is positive-definite, the matrix J + \T\E is invertible for any real 
A. Note that for negative A the function H\ is complex. If we want a real hamiltonian, we 
must take the real part of H\ (while the complex part is a Casimir function). 

The formula for H\ given here is different from the one in [44H46| . The difference is a 
Casimir function. 

3.2 The bad set 

Since we know that our system is bihamiltonian, we can apply the construction of |50| . 
According to this construction, the first thing we should do is to describe the set Bad, 
i.e. the set of points in which the rank of all brackets drops. This is the set where the 
construction of J5Q does not work. 

Proposition 3.3. Let M be a regular relative equilibrium. Then M £ Bad if and only if 
dimKerA/ > 2. 



11 



The proof can be found in Section 4.4 The construction of Section 1.4 of |50| allows us 
to obtain an involutive system of integrals in the neighborhood of any point M £ Bad. But 
in our case there are globally defined integrals (because Casimir functions of all brackets 
of the pencil are globally defined). It is easy to check that the global system of integrals 
locally coincides with the local one (which is not always the case, see Example 2.4 in [50)). 
Therefore, it is possible to apply the theorems of [50] to these global integrals. 

However, we still cannot say anything about the points which belong to the Bad set 
(though our global integrals are defined on this set as well). 

3.3 Complete integrability 

Liouville integrability of Euler- Arnold equations was proved in |20| . This can also be easily 
done using Bolsinov criteria [Hj of integrability of bihamiltonian systems. In both cases 
integrability is established on almost all coadjoint orbits. If we want to prove integrability 
on a given coadjoint orbit, we may apply another argument: 

Proposition 3.4. Let T be an analytic involutive system on a symplectic manifold M 2n . 
Suppose there exists a non-degenerate singular point of T . Then T is complete on M 2n . 

This may sound strange, because non-degeneracy is usually defined for systems which are 
already known to be completely integrable. But actually the definition of non-degeneracy 
works for any involutive system of integrals. And if the condition of non-degeneracy is 
satisfied, then the system of integrals is automatically complete. 



We will give a simple criteria (see Section 4.1l which makes it possible to check non- 
degeneracy of a zero rank singular point. Therefore, we have a simple sufficient condition 
for completeness on a given coadjoint orbit. 

4 Non-degeneracy 

4.1 Non-degeneracy and type theorems 

Recall that a zero rank singular point of an integrable system on a symplectic manifold is a 
point in which the differentials of all integrals vanish. 

Let M £ Bad. Take any bracket P a of the pencil (|8J) such that rank P a (M) is equal to 
the rank of the pencil. Let 0(M, a) be the symplectic leaf of this bracket passing through 
the point M. 

Theorem 6 (Bolsinov, Oshemkov [49]). M £ Bad is a zero-rank singular point for the 
system restricted to 0(M,a) if and only if there exists an orthonormal basis such that J is 
diagonal and M is block- diagonal with two-by-two blocks on the diagonal. 

In other words M £ Bad is a zero-rank singular point if and only if it is a regular 
equilibrium point. 

The following theorem gives a condition for non-degeneracy of a zero rank singular point. 

Theorem 7. A zero-rank singular point M ^ Bad is non- degenerate for the system re- 
stricted to 0(M, a) if and only if the parabolic diagram of M is generic. 

Remark 4.1. Note that a is chosen in such a way that rankP a (M) is equal to the rank 
of the pencil. If this is not the case, it is possible for M to be non-degenerate even if 
the parabolic diagram is not generic. For example consider diagram 7 on Figure [3] The 
corresponding singular point is degenerate on all four-dimensional symplectic leafs passing 
through it. However it is non-degenerate elliptic on the singular two-dimensional symplectic 
leaf. This argument can be used for proving stability of the point. 
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Remark 4.2. A sufficient condition for non-degeneracy is given in [49]. In our terms it 
means the following: 

• The parabolic diagram is generic. 

• For each A there is no more than one intersection point on the parabolic diagram with 
x coordinate equal to A. 

Theorem 8. The Williamson type of a non-degenerate zero-rank singular point M <£ Bad 
is (k e ,kh,kf), where 

• k e is the number of real intersections on the parabolic diagram in the upper half-plane 
plus the number of intersections at infinity, 

• kh is the number of real intersections in the lower half-plane, 

• kf is half the number of complex intersections. 

Scheme of the proof of Theorems [7l [8] 

According to the general scheme of [50] to prove non-degeneracy and find the type of a 
singular point M we should do the following: 

• Find those A, for which the rank of P\(M) drops, i.e. describe the set A(M). This is 
done in Section [4.41 

• Check that the pencil is diagonalizable at M. This is done in Section [4. 5[ 

• Linearize the pencil, check that linearizations are non-degenerate and find their type. 
This is done in Section \4. 61 

• Collect all this together. This is done in Section [4. 7| 



4.2 Examples 

Now we shall discuss examples. 

Example 4.1 (Three-dimensional body). Parabolic diagrams for a 3-dimensional body and 
are illustrated on Figure[2] Rotations around the long and the short axes are elliptic singular 
points (diagrams 1,3). Rotation around the intermediate axis is a hyperbolic singularity 
(diagram 2). 

Example 4.2 (Four-dimensional body). Parabolic diagrams for a 4-dimensional body are 
illustrated on Figure [3] The first case from Example |2.2| corresponds to a center-center 
singular point (see diagrams 1,2). The second case corresponds to a center-saddle singular 
point (see diagrams 3,4). In the third case everything depends on the ratio of angular 
velocities. As we change the angular velocity of rotation in the "inner" plane, the following 
bifurcations occurs: 

Center-center (diagrams 5,6) — >• degenerate (diagram 7) — ► focus-focus (diagram 8) — >• 
degenerate (diagram 9) — > saddle-saddle (diagram 10). 

4.3 Notations 

Let us fix the notations which will be used throughout the proofs. We'll mainly discuss 
regular equilibria (except for the proof of Theorem Eb. Therefore, we may always assume 
that there exists an orthonormal basis such that J is diagonal, while fi and M are block- 
diagonal with two-by-two blocks on the diagonal. 

Let us denote by Xi the diagonal elements of J in this basis. Note that this means 
that Xi are possibly different for different equilibria. However, they are unique up to a 
permutation and coincide with the eigenvalues of J. 
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By Wj's we will denote the non-zero entries of matrix Q, i.e. 



( ° 


Wi 






-Wi 
















Wr, 






- w m 






V 



0/ 



By mi = (A28-1 + )\2i)LOi we wm denote the entries of the matrix M. Mi will stand for 
diagonal two-by-two blocks of M i.e. 



M, 




-rrn 







n will always stand for the dimension of a body, m - for the number of non-zero cj^'s (i.e. 
for the number of two-dimensional planes in the decomposition (|5J) ) . 

For a fixed A let A — J 2 — \E if A 7^ 00 or A — E otherwise. By g^ denote the diagonal 
entries of the matrix A. Clearly Oj = A 2 — A if A 7^ 00 and ai — 1 otherwise. 

It will also be convenient to define the value of Xi( x ) given by H at the point infinity: 

< \ 1 1 



w 2 (A 2 i-i + A2;) 2 m 2 



This agrees with the definition "Two parabolas intersect at infinity, if their leading coeffi- 
cients are equal." 

Pa will stand for the Poisson tensor of the bracket { , }a (given by fflSJ) ) at a given point. 

By sgrad/ we will mean the hamiltonian vector field generated by /: 

sgrad / = Pood/. 
Eij will always denote the matrix with a one in position (i,j) and zeros elsewhere. 

4.4 Description of A(M) 

Let M be a zero-rank singular point. Let us find a basis such that J is diagonal and M is 
block-diagonal. Introduce the following subspaces: 

• K C so(n) is generated by p2i-i,2i — i?2i,2i-i where i — 1, . . . , m and all Eij — Eji for 
2m < i < j < n. 

• Vtj C so(n) is a subspace generated by £21-1, 2j-i — Eaj-i,a%-i, Eu-xjij — E2j,2i-i, 

E2i,2j-1 — E2j-l t 2i, E2i,2j ~ Eaj,2i- 

• Wij C so(n) is a subspace generated by E-2i-i,j — Ej t 2i-i, Eii,j — Ej : 2i. 
We have a vector space decomposition 

so{n) = K® Vn® Wij. 

l<i<j<7n l<i<m, 

2ra<j<?i 

Proposition 4.1. The space K belongs to the common kernel of all brackets of the pencil 
at the point M . All spaces Vij, Wi are pairwise orthogonal with respect to all brackets of the 
pencil at M . 
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Proof. This is a simple straightforward computation. 



□ 



Therefore, the rank of a bracket of the pencil drops if and only if this bracket is degenerate 
on one of Vij or Wij. Let us calculate our brackets on these spaces. 

Identify Vij with the space of two-by-two matrices and Wij with R 2 . Let Mi, . . . M m be 
two-by-two diagonal blocks of M. Let A = J 2 — \E if A 7^ 00 and A — E otherwise. Write 
down A as 

/ A 1 \ 



A 



2n 



where Ai are two-by-two diagonal matrices and en are numbers. 
Proposition 4.2. The form P\ restricted to Vij has the form 

Px(X, Y) = 2Tr (MiXAjY* + M j X i A t Y) 

The form P\ restricted to Wij has the form 

P\(v,w) = -2a,jMi(v,w). 

Proof. This is a straightforward computation. 

Let us now calculate Pa on Va in coordinates. Let 



□ 



A a 



r ij 

a2 S -i 

a 2s 



,X: 



a b 
c d 



,Y: 



e f 
g h 



Explicit calculation shows that 

P\(X,Y) — 2(rrna2j-ic + rrija2i-ib)e + 2{rrijaiid — m i a2j-\a)g J f 
+ 2(rrna2jd — mja2i-ia)f — 2(m,ia2jb + mja2ic)h. 

Consequently X £ Ker P\ if and only if 

iriia2j~ic + rrija2i-\b — 0, 
mja2id — m,ia2j-ia — 0, 
mia2jd — m,ja2i-ia — 0, 
K rma2jb + irija2ic = 0. 

This system can be split into two two-by-two systems and the determinant of both of them 
equals 

det = mj<X2i-ia2i — m,ia,2j-ia,2j. 
Consequently, we have proved the following 
Proposition 4.3. Pa is degenerate on Vij if and only if 



rrija2i-\a,2i — m i a2j-\a2j — 0, 



(10) 



where a s — A 2 — A if A 7^ 00 and a s = 1 if A = 00 . 
The kernel in this case is given by 



X 



amja2i f3mia,2j-i, 
—(3mja2i-i cemiCi,2j-i 



where a and j3 are arbitrary numbers. 
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Now we shall study Pa on Wij. 

Proposition 4.4. P\ is degenerate (and, consequently, zero) on Wij if and only ifX — X^. 

The proof is straightforward. 

Proposition 4.5. The intersection of kernels of all brackets of the pencil is exactly K. For 
almost all brackets the kernel is exactly K . 

Proof. Indeed, only finite number of brackets are degenerate on each Vy and Wij . □ 

Now we are able to describe the Bad set. 



Proof of Proposition \3.3\ Indeed, for almost all brackets the kernel is exactly K, which 
means that all brackets are degenerate if and only if dim A" > [n/2]. This is equivalent to 
the condition dimKerM > 2, q.e.d. □ 

Now we can prove the following 

Proposition 4.6. Let M £ Bad. Then A(M) is the set of x coordinates of the intersections 
on the parabolic diagram of M . 

Proof. P\ is degenerate on Vij if and only if m^a2i-ia2i — rofa2j-i02j = 0. This can be 
rewritten as 

But this means that A is the x coordinate of the intersection point of two parabolas. 

Further, P\ is degenerate on Wij if and only if A = A?. But this means that A is the x 
coordinate of the intersection point of the vertical line with any parabola. □ 

4.5 When is the pencil diagonalizable? 

As a next step, we should check that the pencil is diagonalizable at M. 

Proposition 4.7. The pencil is diagonalizable at a point M £ Bad if and only if any two 
parabolas on the parabolic diagram of M intersect at two different points. 

Proof. Proposition 3.2 of 50. implies that the pencil is diagonalizable if and only if 

C®so{n)/K= Ker (P x |c®„(n)/Jf) • (11) 

AgA(M) 

We have 

€®so{n)/K= C igi Vij© C® Wij, 

l<i<j<m l<i<m, 

Since all the summands of this decomposition are pairwise orthogonal with respect to P\, 



(111 is satisfied if and only if 



C<S Vij = Ker (P A |c®Vy) for l<i< j < m, (12) 

XeA(M) 

C ® Wij = Ker (P A | c ®w y ) for 1 < i < m, 2m < j < n. (13) 

ASA(M) 



Since there is a unique A = Xj such that C <S> Wij = Ker (Pa |c®Wy)i the relation (131 is 
always satisfied. 



Relation ( 12 1 is satisfied if and only if the equation dlOB has two distinct roots, i.e. if 



corresponding two parabolas are not tangent to each other, q.e.d. □ 



1G 



4.6 Linearization 

The last step is to linearize the pencil and to check whether the linearizations are non- 
degenerate. We have 

8t>(n)=K<B Vn® W K 



'13' 



1 < i < j < m 1 < i < m , 

2m<j<n 

The kernel of each P\(M) can be decomposed in the following way 
KerP x =K® % © Wij, 

l<i<j<m l<i<m, 

2m<j<fi 

where Vy C Vy , Wij C Wij . For complex values of A we have Vij C C ® Vij . 

Proposition 4.8. The spaces Vij, Wij are invariant with respect to the adjoint operators 
adx in 0a, where X £ K . 

Proof. Linearization jja in our case is simply the stabilizer of M with respect to the bracket 
[,]a- Consequently, the commutator in Q\ has the form 

[X,Y]x =XAY -YAX, 

where A = J 2 — A_B for finite A and A — E for A = oo. 

Let X € K, Y £ V^. Then it is easy to see that [X, Y]x G Vij, which means that 
adx (Vij) C Vij. But Vj = V^ n KerP A , therefore \4j is invariant. The proof for Wij is the 
same. □ 

Now represent an element X £ K as 

xi 



X 



-»i 



Proposition 4.9. Consider the case when Vij is non-empty. Then the eigenvalues of ad} 
restricted to Vij are ±Vij(X), where 



Vij{X) = V-Xi{XJ(mjXj - rriiXi). (14) 

Proof. This is a straightforward computation. □ 



Remark 4.3. If Vij is non-empty, then XiW — Xj(ty> therefore the formula ( 14 1 is actually 
symmetric. 

Proposition 4.10. Consider the case when Wij is non-empty. Then the eigenvalues of 
adx restricted on Wij are ±/j,i(X), where 

f>i(X) = \/-XiW m i X i- 

Proof. This is a straightforward computation. □ 

Consider the following set R of linear functions on K: 
R = {vij(X) | i,j are such that Vij 7^ 0} U {fii(X) \ i is such that Wij 7^ for some j}. 

Now we see that the following is true 

Proposition 4.11. K is a diagonalizable subalgebra in $\. The set of roots of Q\ with 
respect to K is the set {±^ | £ G R}. 
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Proposition 4.12. The \-lineamzation of the pencil is non- degenerate at M if and only if 
there are no three objects (which may be parabolas or vertical lines) on the parabolic diagram 
which intersect at a point with x coordinate equal to A. 

Proof. To prove our proposition it suffices to show that the linear functions belonging to R 
are independent if and only if there are on three objects on the parabolic diagram which 
intersect at a point with x coordinate equal to A (see Theorem 4 of |50|). 

First note that up to multiplication on non-zero constant the elements of R are 

T>ij(X) — rriiXi — m,jXj, 
Jli(X) — miXi. 

If there are three parabolas intersecting at one point (with x coordinate equal to A), then the 
spaces Vij, Vjk, Vik are non-empty for some different i,j, k. But this implies that R contains 
Vij,Vjk,Vik and since 

Vij + Vjk : — Vik, 

the elements of R are not independent and the linearization is degenerate. 

Now suppose that there is an intersection of two parabolas and one vertical line. Then 
R contains Vij,Jli,Jij and since 

v ij — IM ~ t^j ) 
the elements of R are not independent and the linearization is again degenerate. 

Vice versa, suppose we do not have any triple intersections. Consider some intersection 
point. One of the intersecting objects is necessary a parabola. Therefore R contains either 
some Vij = miXi — mjXj or ju, = niiXi. The coefficient in front of Xi is non-zero in both cases. 
But there are no other elements of R, containing the term irnxi. Indeed, a parabola can't 
intersect any other object at a point with the same x coordinate. Therefore, the element of 
R which is given by each pairwise intersection is independent with the other elements of R, 
which means that R is independent and the linearization is indeed non-degenerate. □ 

Proposition 4.13. Let the X-linearization of the pencil be non-degenerate at M , where A 
is real or infinite. Then the type of Sing(d\H(M)) is (k e ,kh,0), where 

• k e is the number of intersections on the parabolic diagram such that their x coordinate 
is A and their y coordinate is positive; 

• kh is the number of intersections on the parabolic diagram such that their x coordinate 
is A and their y coordinate is negative. 

Proof. Indeed, each intersection of parabolas correspond to the pair of roots ±Vij(X), where 



Vij(X) = y-Xi(XJ(mjXj - miXi). 

The y coordinate of the intersection point is y = XiW an d Vij is real if and only if this 
number is negative. 

Intersection of a parabola with a vertical line corresponds to a pair ±/i;(X), where 



Again, Hi is real if and only if the y coordinate of the intersection given by y — XiW is 
negative. 

We conclude that the number of pairs of real roots equals the number of intersections 
in the lower half-plane, while the number of pairs of imaginary roots equals the number of 
intersections in the upper half-plane. Taking into account Theorem 4 of [50 , this proves our 
proposition. □ 
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4.7 Proof of non-degeneracy and type theorems 

Proof of TheoremFft M is non-degenerate if and only if the pencil is diagonalizable at M 
and all linearizations are non-degenerate (Theorem 5 of [SDj). Pencil is diagonalizable if and 
only if any two parabolas intersect exactly at two points (Proposition|4.7[). All linearizations 



are non-degenerate if and only if there are no multiple intersections (Proposition 4.121. But 
these two conditions together give exactly the condition for a parabolic diagram to be 
generic. □ 

Proof of Theorem [£\ By Theorem 6 of [SO] , the type of a singular point M is k e ,kh,,kf, 
where 

fee = ^2 k e (X), 

AeA(A/)n! 

k h - ^2 k h (X), 

AeA(M)n! 
kf — 2_] kf(\) + - 2_] (dime Ker Pa — corank IT) , 

AGA(M)nK SeA(x), 

ImA>0 

and (k e (X),k h (X),k f (X)) is the type of Sing(d\Tl(M)). 

Let A be real (or infinite). Then k e (X) is the number of intersections with x = X,y > 0, 



fefc(A) is the number of intersections with x — A, y < 0, and fe/(A) = (Proposition 4.131. 
Therefore, 

y^ fee (A) = the number of intersections in the upper half-plane, 

AeA(M)nI 

y^ fcft(A) = the number of intersections in the lower half-plane, 

A6A(M)nI 

AGA(M)nI 

If A is complex, then 

- (dime Ker Pa — corank II) 

is the number of intersections with x — A. Since we count A with Im A > 0, this sum is one 
half of the total number of complex intersections. The theorem is proved. □ 

5 Proof of stability theorems 

5.1 Stability 

Proof of Theorem^ Take a regular symplectic leaf O passing through M. The conditions 
of the theorem imply that the point M has pure elliptic type on O. Therefore, there exists 
an integral / such that 

f(M) = 0, 
d(/|o)(M) = 

and the Hessian of / [o is positive definite at M. 

Since O is regular there exists a coordinate system xi,... Xn in the neighborhood of M 
such that 

O = {xi = 0,i = l,...,k}. 

h 

Now it is easy to see that F — f + Yl x 1 i s a Lyapunov function. D 
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5.2 Non-resonancy 

To prove instability theorems, we need to check that our system is non-resonant, i.e. that 
the trajectories of the system are dense on almost all Liouville tori (see [30] for the precise 
definition of a non-resonant integrable system). This seems to be the most tricky part of 
the whole paper. 

Theorem 9. The system of Euler- Arnold equations is non-resonant, i.e. the trajectories of 
this system are dense on almost all Liouville tori. 

First we will prove several preliminary statements. 

Proposition 5.1. Let M be a regular equilibrium such that the parabolic diagram of M is 
generic and all eigenvalues of M are distinct. Then the eigenvalues of the linearization of 
the Euler-Arnold vector field at M are ±<7y , where 

I" > I x]f + \ 2i -l\2i x\f + \2j-l\ 2 j 

0~;a = 



where 1 < i < j < [n/2] and x^ are two roots of the equation Xi( x ) — Xj(. x )- 
If the dimension is odd, there are also eigenvalues ±Ti, where 

_ 1 /A n + \2i-1^2i - 

where 1 < i < [n/2]. 

Proof. Instead of the linearization of sgradi-/, we can consider the operator D = DhPoo 
(see Section 4 of |5D]). Since P x dH x = P^dH, we have D = D H Poo = D Hx P\- 

Since the parabolic diagram of M is generic, the pencil is diagonalizable at M , and there 
is a decomposition 

C®so{n)/K= CigKer P X (M)/K, (15) 

AeA(M) 

where K is the common kernel of regular brackets of the pencil at M (see Proposition 4.1 
of [50]). 



Since KerP^ is invariant with respect to Dh x P\, the decomposition ( 15 I is invariant 
with respect to the operator D. Now note that 

D I Kcr p x = addir A , 

where ad is the adjoint operator in KerPx- Now, applying Propositions |4.9||4.10| and taking 
into account the formula d9l, we get the formula for the eigenvalues. □ 

Proposition 5.2. Let \2i-1\2i 7^ A2J-1A2J for all i,j. Then o^ and n are linearly 
independent as functions of M. 

Proof. Suppose that 

^ {aljVij + (hjO-lj) + 22 biTi = °- 

Ki<3<[n/2] Ki<[n/2] 

Fix k and choose those members which do not depend on rrik. Their sum does not depend 
on nik, therefore the sum of all other members 

Sk — /2 (^k^lk + a>ik°ik) + 2_j (flkjVkj + a k 3 °kj) + b k r k (16) 

l<i<k k<j<[n/2] 



21) 



does not depend on m k as well. For simplicity denote 

(Tjfc = —akj,a,jk — —dhj 



and rewrite ( 16 1 as 



Sk = 2_^{a ik cr ik + a ik a ik ) + b k r k . 



(17) 



Let m k tend to zero. It is easy to see that 



lim a- 



0, lim r k = 0. 

nik—tO 



Consequently, 



lim S k = 0. 



But S k does not depend on m k which means that S k = 0. 



Now fix I 7^ k. Again, there are two summands in ( 171, which depend on mi. Their sum 



o 11.22 

Jlk — ai k <Ji k + a lk O lk 

must not depend on mi. But this sum tends to 0asmi-> 0, therefore Si k = 0. But 

lim a ik — Ami, lim &i k — Bmi, 



A = ±- 



(A2fc-1 — A2J-l)(A2fc-l — \2l) 
^21-1 + ^21 V (A2fc-1 + A2J-l)(A2fc-l + A2i) 



B = ±- 



(A2fc — A2i-l)(A2fc — \2l) 



^21-1 + ^21 V (A2fc + A2i-l)(A2fc + A2i) 

where the signs depend on the relative positions of A's. Consequently, 

a} k A + af k B = 0. 
On the other hand, 

lim a] k = Cm k , lim af k - Dm k , 

mi— tOO 771;— »00 



(18) 



Consequently, 



C = ±-. 



1 



D = ± 



A2fc-1 + A2fc 

1 

A2fe-1 + A2fc 



(A 2 



A 



)(Aa 



2h-ljy^2l-l 



A 2fe ) 



(Aai-i + A2fe-i)(A2i-i + A2fe) 

(A2i — A2fe_l)(A2( — A2fc) 
(A2i + A2fe_l)(A2! + A2fc) 



a ik C + a tk D — 0. 



(19) 



It is easy to see that AD — BC — if and only if A 2i-i A2i = A2fc-iA2fc. But we assumed 
that this equality is not satisfied. Therefore ( 18 \ and ( 19 \ imply that a] k — af k = 0. Since k 



and I were arbitrary, all coefficients a i '■ vanish. But this implies that b; vanish as well and 



our functions are linearly independent, q.e.d. 



□ 
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Remark 5.1. It is easy to see that if A2i-iA2i = A2J-1A2J for some i,j, then o~\j — crjj, 
which means that the eigenvalues are not linearly independent. 

Now the proof of Theorem [9] follows from the following simple lemma: 
Lemma 5.1. Let /i,.../ n be continuous functions on a manifold N. Suppose that these 
functions are linearly independent on any open subset V C N. Then their values are inde- 
pendent over Z almost everywhere on N . 

Proof of Theorem^ Let us consider an equilibrium M such that Ai < A2 < • • • and the 
parabolic diagram is generic. In the neighborhood of M the eigenvalues of sgrad 1/ are 
linearly independent. Then, by the previous lemma, they are independent over Z almost 
everywhere. Take a point Mi such that they are independent. We claim that our system is 
non-resonant in the neighborhood of Mi. Indeed, Mi is an elliptic point and the Liouville 
foliation on any regular symplectic leaf passing through Mi is locally given by the functions 

Si = Pi + qi, ■ ■ ■ , Sm =Pm + Qm 

in some symplectic coordinates p, q (Eliasson theorem [28]). 

The functions st are action variables and the rotation numbers are given by 

Ci(si,...,s m j- 9fl _. 

It is easy to see that Ci(0) are exactly eigenfrequencies of sgrad .ff. They are independent 
over Z. But this implies that a are independent on almost all tori in the neighborhood of 
Mi and our system is non-resonant. By analyticity, it is non-resonant everywhere. □ 

5.3 Instability 

To prove the instability theorem, we will need the following simple 

Lemma 5.2 (About unstable cone). Let a vector field v on a manifold N vanish at point 
xo . Suppose that the linearization of v at xo has an eigenvalue with a positive real part. 
Then there is an open subset K <Z N such that 

1. There exists 8 > such that all trajectories starting at K leave Us(xq). 

2. The intersection of K with any open neighborhood of xo has non-empty interior. 

Corollary 5.1. Let v — sgradH be a non-resonant integrable hamiltonian system and x be 
an equilibrium point of it. Suppose that there exists an integral f such that sgrad /(a;) = 
and the linearization of sgrad f at x has an eigenvalue with non-zero real part. Then x is 
an unstable equilibrium for sgrad H . 

Proof. Since the linearization of sgrad/ at x has an eigenvalue with non-zero real part, it 
has an eigenvalue with positive real part. Therefore, we can find a set K from the lemma. 
For any e the intersection U E (x) Pi K has non-empty interior. Therefore, we can find a 
non-resonant torus passing through U e {x) n K. The trajectory of sgrad/ lies on this torus, 
therefore this torus will leave Us(x). But since the torus is non-resonant, all trajectories of 
sgrad H are dense on it and will leave Us (x) as well. Therefore, x is an unstable equilibrium, 
q.e.d. □ 

Proof of TheoremlM By Proposition 4.4 of [50] the following sets of operators are equal 

{DfPoo !KcrpJ /6Jid/eKcr p m = {ad 4 } ? gK, 

where adg is the adjoint operator in g>,- 

Conditions of the theorem imply that for some A there is £ G K such that the operator 
ad^ has an eigenvalue with non-zero real part (see the formulas for the eigenvalues given by 
Propositions |4. 9| |4. 10[ ) . Therefore, there exists an integral / such that DfPao also has such 
an eigenvalue. But DfPoo is dual to the linearization of sgrad/. Now it suffices to apply 
Corollary |5T| □ 
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5.4 Instability of exotic equilibria 

Now we shall prove that all exotic equilibria are unstable. 

Proposition 5.3. Let M be an exotic equilibrium. Then there is A such that sgrad .Ha 7^ 0, 
where 

H x = -^{(J + ^/\E)- 1 Q{J + V\E)- 1 ,M). 

Proof. The vector field sgrad H\ has the form 

til = [{J + vE)- 1 Q.(J + vE)- 1 ,M], where v = \/A. 
Suppose that M is an equilibrium point for sgrad H\ for all A. Then 

= [(j + iyE)- 1 n(j + uEy 1 ,jn + nj] = (,/ + uEy 1 ^ 2 - o. 2 (j + vE)~ x + 

+ ( J + vEy l n(J + i/_B) _1 fi(J - vE) - (J - vE)Q.(J + vE)~ l Q.{J + vE)' 1 . 

Since M is an equilibrium point of the body, 17 2 commutes with J. Therefore, it also 
commutes with (J + vE)~ J . Consequently, 

(J + vEy x Q.{J + uEy^iJ - vE) = (J - vE)Q,{J + vE)~ x n{J + vE)' 1 . 

But this equality means that the matrix 

(J + vE)' 1 ^ + vEy l Q.(J - vE) 

is symmetric. Let us denote by u>ij the entries of the matrix Q. \t are, as usual, diagonal 
entries of J. Then the symmetry condition can be written as: 



Afc — V -s--\ LOijLdjk Aj — f ^— v (jJijUJjk 

\i + v Is Aj ■ + v ~ A fc + v ^r 1 Aj • + v 

3 3 

Since all eigenvalues of J are distinct, this implies 



Y ^± = for all i + k. 



A, + v 

3 

This equation should be satisfied for all v. But this means that 

LOijUljk — 

for all j and all distinct i,k. Therefore, we can bring SI to the block-diagonal form with 
two-by-two blocks on the diagonal by permuting basis vectors. Such a permutation will 
preserve diagonal form of J. Consequently, M is not exotic, but regular. Contradiction. □ 

Proof of TheoremUl By the previous Proposition we can find an integral / such that 

sgrad /(M) /0 

for a given exotic equilibrium M. Obviously, the trajectories of sgrad/ leave sufficiently 
small neighborhood of M. Therefore, Liouville tori leave this neighborhood as well. Since 
our system is non-resonant, it's trajectories are dense on most Liouville tori and will also 
leave the neighborhood, q.e.d. □ 



Remark 5.2. We used Proposition |5.3| to show that exotic equilibria are not rank zero 
singular points. For equilibria not belonging to the set Bad this follows automatically from 
Theorem [6] 

Remark 5.3. Since sgrad H\(M) 7^ for any exotic equilibrium M and some value A, 
exotic equilibria are not isolated on the symplectic leaves of the so(n)-bracket, but form 
smooth families. The dimension of such a family essentially depends on the sizes of the 
blocks Ai, entering formula Q. 
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